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ABSTRACT 



Computations indicate that when utilized in an idealized 
ramjet at a flight Nach number of 2.5 and an altitude of 
37,500 feet, a constant area axial flow nuclear heat exchanger 
will be expected to produce 9,500 pounds of thrust with an air 
inlet Kach number of 0.2, a weighted radial temperature aver- 
age of 2000°? and a frontal area of 7.7 square feet. Molyb- 
denum alloy with a diailicide coating was selected as the 
container material. Sandwiched between the molybdenum was a 
fissionable alloy of highly enriched uranium and zirconium. 
Prom a nuclear aspect the unreflected heat exchanger becomes 
critical at a diameter of ipO inches and a length of 60 Inches 
with a total weight of 8,052 pounds which includes 3,622 
pound 3 of 90 per cent enriched uranium. 
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introduction 

That a constant-area axial-flow compact nuclear heat 
exchanger as envisioned in Figure 1 when installed in an 
idealized ramjet woxild produce significant thrust was the 
basis of this thesis. Feasibility conclusions were limited 
to the information revealed by heat transfer and nuclear 
analysis under one assumed steady state condition. The 
problems of nuclear control, shielding, thermal stresses, 
and performance under different inlet Mach numbers wero not 
considered. 

The nuclear heat exchanger as shown in Figure 2 incor- 
porates thin concentric cylinders which are the solid fis- 
sionable containers. The alloyed fissionable phase is 
sandwiched between two layers of non-f issionable structural 
material so that the enriched uranium fuel cores are sur- 
rounded by a continuous network of structural metal for 
strength and containment at high temperatures. Air flows 
within the annuli and acts as the coolant and thrust agent. 

It is recommended that the nuclear flux be flattened by a 
central control core and a peripheral reflector; however, the 
beneficial effects of a reflector and central control ooro 
were not considered. 

Prom the assumed inlet condition, various maximum cyl- 
inder or wall temperatures were selected arbitrarily for 
thru3t results. That the high temperature fuel elements and 
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structural materials wore possible install urgically wer c next 
Investigated. Whether or not an unreflected homogeneous 
reactor of this conf iguration can be made critical was deter- 
mined by using one group diffusion theory. 
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LIST OF SYMBOLS 



w 



H 

k 

K 



Cross section area normal to flow, (ft^), or the rsaas 

number of an element (gny'gm atonic weight)* 

Wetted area, (ft 2 ) . 

- 2 , 



Buckling, (cm ). 

Local velocity of sound, (ft/sec). 

Specific heat at constant pressure, (BT0/lb °P). 
Diffusion coefficient for flux, (cm). 

Hydraulio diameter, Dg - D^, (ft). 

Inside diameter, (ft). 

Dg Outside diameter, (ft). 

2 3 

-W / Neutron leakage, ( n / exsr sec). 

f Friction factor defined as 2g -r v /? v -- 

p 

g Acceleration of gravity (3 2.2 ft/sec ) . 

p 

0 Mass velocity (lb/sec f t ) . 

h Coefficient of heat transfer between fluid and surface. 



(BTU/hr ft 2 *P). 
h ft based on (T 0 w “ T 



o,aw 



). 



h r baaed on logarithm mean of (T - T) . 

u Q|W 



h based on (T . - T) . 



at 

h based on (T 
s 



o,w 



o,w 



- V- 

Extrapolated distance, (cm). 

Neutron multiplication factor. 

Thermal conductivity, (BTU/hr ft 2 •F/ft). 
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L Length of heat exchanger heating element In the 

longitudinal direction, (ft), 
ra Haas rate of flow, (lb/a«o). 

M Mach number or molecular weight (gas/gm aolecular 

weight) • 

n ?fumber of neutrons per cir*. 

N Humber of nuclei por cra^. 

O') 

Avogadros number ( 6.02 x 10 ^ nuclei/gra atomic weight). 
"pR Prondtl number, C pt i/K. 

Kre Reynolds number, 9 vD fl /p. 

Stanton number, h^/C^G. 

P Absolute stream pressure (lb/in ). 

Pq Isentropic stagnation pressure (Ib/in ). 

q Heat-transfer rate (BTU/hr). 

Q Volumetric heat source which la equal to E 

(BTO/aec cirr^). 
r Radius (ft). 

T Absolute stream temperature (°R). 

Thrust (lb). 

Tq Stagnation temperature ( # R). 

T 0 aw Adiabatic stagnation wall temperature (*R). 

T q m Stagnation wall temperature (°P) , 

▼ Velocity (f t/aeo) . 

V Volume (ft^ or c*3). 

W Weight (lb). 

x Distance along heat exohanger heating elements (ft). 
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Subscript a 
b Bulk (T) * 

f film t f = (t^ + t w )/2. 

w Wall. 

0,1,2 Station rofcrenco3 
Greek 

a Definition of vJi/tf 

if Ratio of specific heats of gas, C p/ C v ** 

Ah Change in enthalpy, (BTU/lb). 

Tj Number of neutrons produced per absorption in the fia- 
slonablo isotope. 

/( Definition of 

Transport mean free path, (ca), 

Absolute viscosity of the fluid, (lb/seo ft). 

P o Average cosine of th© scattering angle per neutron 

collision in the laboratory system, 
r' Number of neutrons produced per fission in the fission- 
able isotope. 
p Density, (Ib/ft 3 ). 

p 

<7“ Microscopic cross section (cm ). 

2 

<T microscopic absorption cross section (cm ). 

M 

2 

<r # microscopic scattering cross section (era ). 

<T microscopic fission cross section (cm ). 
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S 



T w 

0 

72 



Macroscopic cross section (cm* 1 ), (H(T). 

S a macroscopic absorption cross section (cm* 1 ), 
g macroscopic scattering cross section (cm* 1 ), 
macroscopic fission cross section (cm* 1 ). 
Number of neutron absorptions / car* sec. 

Tractive force at wall due to fluid friction (lb/ft^). 
Flux in neutrons / cm sec. 

Implacian operator. 
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LIMITING CONDITIONS 

To obtain a realistic operating condition for tho air- 
cooled nuclear heat exchanger the following conditions were 
assumed: 

M 0 * 2.5, h * 37,500 ft, T * 392®H, P - 3*09 lb/in 2 

Godsey and Young (1, p.221) stated that the total temperature 
through the inlet diffuser depends only on flight Mach num- 
ber and ambient temperature and is a constant regardless of 
diffuser losses if the process is considered adiabatic. 
Therefore the relation 

T 0,0 - T 0.1 = t 0.2 » TU + «o Z ) u> 

nay be used where the station numbers correspond to those of 

Figure 1* Godsey and Young (1, p.221) stated that the total 

pressure at the heat exchanger inlet, P n „, was dependent upon 

u,<: 

the flight Mach number, ambient pressure, free stream diffu- 
sion losses, duct friction losses, end duct diffusion losses. 
Tho last three effects may be grouped under the terra ram ef- 
fectiveness or pressure recovery effectiveness. From data 
presented by Godsey and Young (1, p*86), the ram effectiveness 
at Mq * 2.5 was given as approximately 0.98; however, subsonic 
diffusion was assumed to be loss free. For a more realistic 
condition a value of 0.90 was chosen. The above ratios may be 
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written out in the following expression for stagnation pres- 
sure at the heat exchanger inlet: 



0,2 



p (1 + r- 1 M z ) ’A* 1 x P °.2 actual 
P 0 U 2 M 0 ’ P 0>2 ldeal 



(la) 



Solutions to Equations 1 and la give ^ - 8314.°^ and 
p 0 a 4.6,5 lb per sq in, for inlot conditions to the heat 
exchanger . 

Prior to converting the above inlet stagnation pressure 
to stream pressure, it is necessary to examine the limiting 
conditions for heat exchange. Assuming a uniform cross 
sectional heat exchanger area under a given inlet stagnation 
temperature, which in thic case is 881f*R, and various inlet 
Waeh numbers, it is desired to exchange heat to the extent 
that the stagnation temperature at the exit of the heat ex- 
changer be raised to a maximum value. By use of methods out- 
lined in Reference ( 2 , p. 1^.210 the results of various inlet 
Mach numbers are shown on Figure 3 which reveal that when a 
oertain maximum temperature occurs, thermal choking oven 
under frictionlese conditions exists in the heating process 
and definitely limits the capacity of the constant area 
combustor or heat exchanger. Stated another way, if at any 
position downstream for a constant area duct, the stagnation 
temperature is equal to Tq ^ maximum, the flow is choked 
since tlie Kach number is unity. Any addition of friction 
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will reduce the downstream interval for ?q ^ me-ximum to occur. 
With M 2 ■ 0.2 and T Q -j maximum achieved, 5l00 o R, the stagna- 
tion pressure would suffer a 20 per cent loss due entirely to 
the heating effect. 

Without becoming involved with a specific shock cone and 
diffuser design, the inlet ^ach number was set at 0.2, which 
resulted in the following stream conditions: 

P 2 = k$.l lb/in 2 
9, « .1385 lb/ft 3 
v 2 * 291 ft/sec 
T 2 * 08O°R 

G * i}.0.4 lb/ft 2 sec 

As shown by Figure 2, the heat exchanger elements are 
thin hollow constant concentric area right cylinders whioh 
are adaptable to one-dimensional compressible gas flow treat- 
ment without introducing serious errors. This assumes that 
in any one annulus the change of stream properties in the 
direction of flow is much larger than in the radial direction. 
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HEAT TRANSFER ANALYSIS 

HcAdama (3# p.289) correlated the experimental data of 
Kays and London Ut) and baeause compactness of an air-born© 
heat exchanger ia an important parameter, a comparison of 
heat transfer coefficients per unit volume for a given fric- 
tion power per unit volume showed ruffled fins and louvered 
plate fins as being the most effective. It appears that the 
highest ratio of heat transfer to pumping power is obtained 
when form drag ia absent, or when all the drag is due to skin 
friction which is the case in common for flat plates or fins. 
Even though externally finned tubes are of value in increas- 
ing the rate of heat transfer, friction losses become pro- 
hibitive for this typo of installation. For those reasons 
thin concentric hollow cylinders were chosen as the fuel ele- 
ments of tho heat exchanger and were restricted to constant 
longitudinal cross section area for nuolear simplif lcation. 

Development 

If heat loss and longitudinal thermal conductivity are 
neglected, the rise in the enthalpy of the air per second is 
related to the effective coefficient of heat transfer from 
an annulus by a differential equation which is developed In 
Appendix C as 
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t (D 2 2 - O^lfT O p dT 0 = h e *(0 2 + D^Tq^ - T 0fBM ) d*. (2) 

The particular temperature differential is shown schematically 

in Figure 4* 'The temperature increment for heat transfer is 

the difference between the actual wall temperature, M , and 

the temperature T n ... which the wall would assume for the 

u # aw* 

same flow conditions but with zero heat flux. McAdams et al. 
(5) found that the effective heat transfer coefficient, h ft , 
expressed in terms of the Stanton number, N^, was independent 
of the temperature potential, Tq w - Tq, when based on 
T q w - T o,aw whil * thc coefficients h^ and h # were not inde- 
pendent. 

The Reynolds analogy with lf pR * 1 relates the Stanton 
number to the friction factor by 

h . 

W 

From experimental data for a moderate AT, — 60, and 

5,000<N R <200,000, the following equations from McAdams 
(3, pp. 1^5,219) my be combined 



f « •046/(DG/p f )* 2 
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(DG/p f )* 2 



(4) 



( 5 ) 










f •,**!«*- 












Ml 









mi 




. 



12 



to produce 





Af 



~T 



(6) 



where, by definition, ^ • Since Equation* 3 and 6 
ere essentially the same, the limitation to & moderate tempera- 
ture change my be removed without introducing serious errors 
when using the previously defined temperature potential. The 
following relationship, where specific heat is considered a 
constant, is used for the remainder of the analysis; 



h 



o 



c p° 




(7) 



The friction factor of Equation 7 is set by Equation 4 which 
is dependent upon Reynolds number and the degree of roughness. 
Since mass velocity, 0, and specific heat are essentially 
constant, the friction factor may be considered to depend upon 
‘ hydraulic diameter, viscosity, and degree of roughness. An 
average value of *f * *008 was obtained from Equation 4 for 
1/2 M < D e < 1-1/2 M and 140Q°B < <1600°R. The variance of the 

friction factor within the above range of hydraulic diameter 
may be compensated for by added degrees of roughness j however. 
Cope (6) found that with water in the turbulent region even 
though friction was six times that for smooth pipes, the heat 
transfer was only 20 to 100 per cent greater than for smooth 
tubes. McAdams (3, p.223) states that because of a very 
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small temperature change the heat transfer coefficients ob- 
tained were somewhat uncertain# Correlation with air may not 
apply due to different Prandtl number effect j however, there 
is a possibility that with small variation in hydraulic diam- 
eter the friction factor could remain a constant by varying 
the degrees of roughness# 

Experimental data have shown that for turbulent flow in 
rectangular ducts and annular spaces friction may be evaluated 
from data for circular pipes using a hydraulic diameter shown 
in Figure 2 as 



A combination of Equations 2 and 7 with Equation 8 results in 



If a gas stream of uniform temperature is brought to rest 
lsontropically, as at the stagnation of a blunt body, the tem- 
perature rise for an ideal gas will be 




(8) 




(9) 
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a T(1 ♦ -S—i M 2 ) 



( 10 ) 



For convenience a temperature recovery factor as shown in 
Figure la defined as 
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T - T 
u.aw 

— p 

l 0 - T 



m _ T 

^Q.aw 

V /2g JC 



(ID 



By & combination of Equation 10 and 11, the stream temperature, 
T, may be eliminated to produce 



T 

Q.aw 



1 



(1 “ 

1 ♦ 



N RF } 4 ^ ^ 

23 "? 



( 12 ) 



Subatitution of Equation 12 into 9# yields 



dT 0 



ft 



- l ♦ 



(l - 



Wpp) 



-t 1 M 2 



1 + 



f -1 **2 

a n 




( 13 ) 



The generalized differential Mach number equation of 
Reference (2, p*4»12) eimplifiea to 



^^2 (1 4- M 2 ) dT Q 



■f 



M 2 (l + M 2 ) 

1 - M 2 




( 14 ) 



with negligible structure drag forcos and constant specific 
heat, molecular weight, area and mass flow rate* 
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Equation 13 was substituted Into Equation lij. which gave 

M 2 ) 



dM 2 « 2f 



dx 

D_ 



( 




< 








< (i + d m 2 )A 


o.w , t a BF 2 

_ X ♦ JLZk 

T i » 2 n 


♦ 2r m 2 



(15) 



For a turbulent boundary layer over a flat plate, 

McAdams (3, p.311) recommends that the temperature recovery 
faotor and the Frandtl number be related by 

®RF - »PR V3 • < 16) 

McAdams et al. ($) reported recovery factors for subsonic flow 
of air near the outlet and of a long smooth tube wore in good 
agreement with those predicted by Equation 16. Appendix A 
shows the procedure for solving Equations 13 and 1$ for one 
Mach number and stagnation temperature at an assumed position 
downstream. From these calculations it was evident that the 
contribution of the recovery factor term was negligible 

for Mach numbers up to approximately 0*8j however, the re- 
sulting stagnation temperatures will be nearly l,$ per oent 
higher than those which have incorporated this term. 

For the remainder of this analysis the recovery factor 
and Prandtl number were assumed unity which reduced Equation 
13 to 
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dT 



0 




(17) 



It Is of Interest to note that Equation 17 could be expected 
from experimental data by McAdams (3, p.3l4) because the heat 
transfer coefficient, h , decreases and approaches the constant 
value of h for temperature differences greater than 100°F. 

For a high temperature potential, h approaches h and may be 
used directly in Equation 1. 

Equations 14 and 17 were integrated approximately over 
short intervals of Mach numbers by using the procedure out- 
lined in Fef orence (2, 4*60) • Appendix B contains the step- 
wise procedure used and Tables 1-3, Appendix B, have the com- 
plete data for stagnation and stream properties as a function 
of L/D ratios for various longitudinally constant wall tem- 
perafcures. Stagnation torperatures and pressures versus L/D 

© 

ratios are shown in Figure 5>* 

However, from a nuclear standpoint tho longitudinal flux 
of fuel element cylinders will follow a sine distribution with 
the maximum flux produced at the midpoint. This does not 
imply that the maximum wall temperature will occur at the 
longitudinal midpoint. By use of 



T 0,3 " T 0,2 




1 - sec a max 



(18) 
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which was developed In Appendix C and by using the previously 
selected temperatures for maximum wall temperatures, the 
stagnation temperatures versus lyt^ ratloB were computed as 
shown in Table 4. of Appendix C and plotted in Figure 5» Stag- 
nation prossuros for variable wall temperatures were obtained 
by interpolation from previously calculated pressures. 

An laentroplc expansion using air data from Keenan and 
Kaye (7) was assumed In order to obtain the change of enthalpy 
from stage 3 to !* of Figure 1. The Brayton cycle as shown in 
Figure 6 represents an overall view; however, at various L/D 

© 

ratios, the change of enthalpy was computed using an ambient 
pressure of 3*09 lb per sq in* Figure 7 shows the variance 
of enthalpy for both constant and variable wall temperatures* 
The right hand ordinate of Figure 7 related change of enthalpy 
to pounds of thrust per aq ft by 




where v Q * 2.5c and v^ * 223. * 

Discussion 

The calculated values of stagnation temperatures and 
pressure as shown in Figure 5 do not show tho contributing 
effects on heat transfer and friction caused by the length 




1 ? ♦ 1 



IQ 



to diameter ratio and the effects of sudden enlargement or 
contraction of flow. The former effect la small since the 
L/D ratio for this heat exchanger should be greater than 60; 
the latter effect may still be appreciable even though the 
peripheral thickness of the cylinders are in the order of 0.05 
to 0.07 in. The specific heat ratio was set at 1*4# and the 
friction factor, specific heat, and mass velocity were con- 
stant throughout the beating process. The cosine flux distri- 
bution in the longitudinal direction Is justifiable for a 
finite cylinder. 

The temperatures obtained by normal forced convection 
methods for an assumed Tq w * 246 0°R are shown In Figure 5> for 
comparison. The temperatures show good correlation with stag- 
nation temperatures obtained by the Reynolds analogy but are 
considerably higher than the stream temperatures. Tho forced 
convection formulas will not ahow any decreasing trend in 
stream temperatures or a maximum stagnation temperature as 
heat Is continually added as will the Reynolds method for com- 
pressible fluid which will eventually produce choked flow at 
a certain L/C) ratio downstream. 

Figure 7 shows that the effect on thrust reduction was 
quite severe for a 2000*F constant wall temperature versus a 
2000°F maximum wall temperature ; the thrust per square foot 
was reduced by 18 per cent. For a uniform velocity at the 
exit of stage 4 # a constant change In enthalpy is required. 

If the neutron flux was such as to produce a series of maximum 
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wall temperatures from 1600 to 220;. °F and with a constant 
enthalpy change of 200 BTtf per lb imposed, the L/D^ ratio 
would vary from 47 to 100. Since all cylinders ore desired 
to have the sane length for nuclear considerations, the hy- 
draulic di&nofcera would vary by 47 per cent with decreasing 
hydraulic diameters as radial distance is increased. This 
diameter variance would still be subject to the limitation 
Imposed on the hydraulic diameter due to a constant friction 
factor. 

To increase the thrust which Is based on the heat ex- 
changer cross-sectional area, each fuel cylinder should oper- 
ate at the maximum enthalpy change as shown by Figure 7 and 
the thrust would be a weighted average depending on the 
frontal area of the heat exchanger. With a constant length 
assumed the hydraulic diameter would very by only 12 per cent. 
However, each annulus would be operating at a different veloc- 
ity level which may create internal vibration within the ram- 
jet pod. With a weighted radial average of 2000*1? for the 
cylinder temperatures encompassing 7.7 sq ft of heat exchanger 
cross-sectional area, from Figure 7 the resultant thrust would 
b© approximately 9,500 lb at a flight Mach number of 2*5 at 
37,500 ft. 

With inlet Mach numbers greater than 0.2, the L/D ratio 

© 

would be substantially reduced. Prom Reference (2, p.431). 
with M 2 * 0.5, f » 0.005, T Q 2 * 1000°R, P Q ^ 2 * 10 lb/in 2 , and 
^0,w * 2000°R and l500°R, choked flow occurred at an B/D^ ratio 
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of 20 and 35, respectively. To resist the pancaking effect 
which may reach the limit for nuclear criticality, the inlet 
Mach number appears to have a definite upper limit when heat- 
ing la desired. 



Conclusion 

For the assumed inlet conditions, cylindrical fuel ele- 
ments and wall temperatures, the amount of hoat transferred 
was substantial enough to produoe 1,000 to 1,400 pounds of 
thrust per sq ft if heat exchanger frontal area with a maximum 
wall temperatures range from 1600 to 2200*P at a flight Kach 
number of 2*5 and an altitude of 37,500 ft. With a weighted 
radial average of 2000°F for the cylinder temperatures encom- 
passing 7.7 sq ft of heat exchanger area the resultant thrust 
was 9,500 lb. 
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FUEL ELEVEN? CONSIDER AT I OH 

Due to the absence of effective moderation within this 
nuclear heat exchanger, the unit may be considered as a fast 
nautron reactor and as such the restrictions normally imposed 
on structural materials by high thermal neutron absorption 
may be removed. Whether any structural material now In exist- 
ence can bo used to contain, protect, and support the fission- 
able phase for the high temperatures required ia the subject 
of the following development and discussion. 

Development 

The nuclear heat exchanger as shown in Figure 2 incor- 
porates thin concentric cylinders which are the solid fission- 
able containers. The enriched uranium fuel cores are to re- 
main a solid to prevent serious containment problems. The 
thin sheets of the fissionable phase are sandwiched between 
the inner and outer layers of the non-fisslonable structural 
metal which predominates in volume so that the alloyed en- 
riched uranium fuel sheets are surrounded by a continuous 
network of structural metal for strength at high temperatures. 

The following are some criteria that were considered for 
the selection of fuel elements and structural material: 



* 
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A. Fuel element criteria 

1. Melting point above 2S00®P 

2. Heutron's energy maintained 

3. Mechanical properties similar to the container 
material during fabrication and high tempera- 
ture operation. 

4* High thermal conductivity and low specific 
heat for good heat transfer 

£>. Ease of fabrication 

B. Container material criteria 

1. High temperature mechanical properties under 
static and small radial dynamic loads 

2. Resistance to oxidation at high temperatures 

3. High resistance to thermal shock 

4.. Koutron energy maintained 

£. Ease of fabrication, rolling, and weldability. 

Discussion 

The melting point of natural uranium is 2072*?; for the 
required temperature range of 2000 to 2200*?, the unalloyed 
uranium would become a liquid and present serious containment 
problems. In order to increase the molting point, the various 
binary diagrams of uranium alloys of Sailer (8, p. 413-414) 
revealed that a uranium-zirconium system of 20.7 per cent zir- 
conium by atomic weight would have a melting point of 245>0°F. 
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Elements with less than a mass number of I 4 .O were not consid- 
ered because of the Increased slowing down effect on the 
neutrons. 

Zirconium alloys at high ter^peraturea suffer a loss of 
strength and a loss of corrosion resistance. As reported by 
Dayton <9, p. U.60) , the strongest alloy tested at 930°F creeps 
at a rate of lCf ^ per cent per hour with a stress of only 
I 4 .OOO lb per sq in. Since a thin sheet is sandwiched between 
two layers of structural material, the loss of strength and 
corrosion resistance are secondary to the requirement that the 
fissionable material remain a solid. The coefficient of 
linear thermal expansion for 20.7 per cent zirconium by weight 
was not available; however, Dayton (9, p#46U3 reported a co- 
efficient for pure zirconium at 932*P for the "a" axis of 
£ 

5.8 x 10 * while that of pure uranium was anisotropic as shown 
by Sailer ( 8 , p.389)* The coefficient of thermal expansion 
for this particular alloy was not found; however, from Sailer 
( 8 , p.234) the composition for 20.7 per cent zirconium by 
weight indicated 3 phase changes for the temperature range up 
to 2200 °P and whether the aniaotropio effect of uranium would 
result in volume changes so sever© as to make bonding or 
cladding impractical was not determined. 

The low thermal conductivities Of both zirconium and 
uranium are partially compensated for by low specific hoats 
and the thin sheet effect of high area to volume ratio. 
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Uranium can be hot rolled without danger of cracking. 
Whether this particular alloy can be rolled into sheets with 
thicknesses as small as 0.03 inch was not determined. 

The container material for temperatures up to 2200®F 
presents a considerable problem when restricted to present day 
material availability. The wrought cobalt-base alloys, such 
as Haynes Stellite-25» Allegheny- Ludlum Steel Corporation 
3-590 and S-816, and Weatlnghouse Electric Corporation Re- 
fract aloy 26 and 70, appear to have excellent corrosion re- 
sistance up to temperatures as high as 2000®P, Eichen and 
Jackson (10, p.556)> in an oxidizing atmosphere. Eowover, 
Kichen and Jackson (10, p.539) reported the ultimate strength 
of Haynes Alloy 25 at 2200®F to be 4,000 lb per sq in- for a 
sheet 0.042 inch in thickness and annealed at 2200®?. Another 
major drawback may be the high residual activity in the fora 
of gamma radiation which is characteristic of cobalt-baae al- 
loys. 

Molybdenum is of interest as a hlgh-tcmperature material 
because of its high aelting point, 4?50 e P, high strength and 
stiffness at elevated temperatures. Northcott (11, p»106) 
reported creep-rupture strength for wrought arc cast molyb- 
denum, .45 per cent titanium, of 25,000 lb per aq in at 2000®F 
for 100 hours. The thermal expansion coefficient as reported 

by Northcott (11, p.26) ranged from 5«1 x 10*^ at 932®F to 

—6 

7.2 x 10 at 3632®P. Thus the thermal expansion of both pure 
zirconium and molybdenum alloys appear to be very similar. 
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However, alloying with uranium might alter the thermal ex- 
pansion characteristic considerably. 

Uncoated molybdenum is useless as a material of construc- 
tion in oxidizing atmospheres abovo 10D0 e F. Norden (12, 
pp. 293-304) reported that so far no molybdenum-base alloys 
have been developed with satisfactory oxidation resistance; 
however, Hioster et al. (13) reviewed progress made with al- 
loying techniques. Sherwood (14, p.26) stated that molybdenum 
specimens protected by coatings 0.003 in thick and consisting 
primarily of molybdenum disilioide have successfully resisted 
oxidation at 3092°P for more than 1000 hours. However, the 
coating must be pore-free in order to prevent serious loss of 
material. Effects of the disiliclde coating on heat transfer 
and t hemal cycling were not included. 

The high thermal conductivity of molybdenum, .32 cal per 
eec cm °C, and low specific heat, .061 cal per gra *C, permit 
it to be heated or cooled rapidly with low resultant stress 
which give it excellent resistance to thermal shock. The heat 
transfer characteristics are excellent due to the high thormal 
conductivity. The neutron moderation effeots are small since 
it has a mass number of 96. 

Aa brought out by Horden (12, p.298) ft chief advantage 
of molybdenum over ceramics, cermets, and some super alloys is 
the variety of forms and sizes available on a production basis. 
By the arc-cast methods, sheets from 0.005 to 0.020 in. thick 
and 36 in. wide are available. Molybdenum has been joined 
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successfully by spot welding, arc welding, soldering, and 
brazing but precis© controls must be used since excess heat- 
ing will result In extremely brittle welds. 

Conclusion 

Molybdenum alloy, 0.45 per cent titanium, with a protec- 
tive dislllclde coating appeared to have the best probability 
of withstanding temperatures up to 2200®P la an oxidizing at- 
mosphere besides being able to withstand thermal shock and 
offer near optimum heat-transfer characteristics. The inner 
fissionable phase of enriched uranium alloyed with 20«? per 
cent by weight of zirconium, even though soft and pliable at 
this temperature range, which could be an advantage as far as 
the container material is concerned, should expel the heat 
produced Internally due to the high surface to volume ratio. 
Whether or not the thermal expansion coefficient of the 
uranium- zirconium would be compatible with molybdenum and that 
a metallurgical bond would be formed with molybdenum was not 
ascertained. 
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NUCLEAR ANALYSIS 

The actual heterogeneous reactor conf Iguration as shown 
in Figure 2 was treated as a homogeneous unref looted system 
using one-group diffusion theory to obtain the critical mass 
of uranium. The flattening of flux and reduction of critical 
mass by use of a reflector were not considered. The neutron 
flux necessary to produce a maximum fuel element temperature 
was calculated. The control mechanism necessary to maintain 
the flux level was not considered. 

Development 

The diffusion equation for thermal neutrons for steady 
state conditions as developed by Glassfeone and Edlund (15, 
p.192) is 

DV 2 ^ (r) - E & fS (r) + pq (F, T th ) - 0 . (20) 

For an Idealized fast reactor composed of atoms of material 
so heavy that neutrons upon colliding rebound without losing 
an appreciable amount of energy, the slowing down density term, 
P<1 (r, T th ), may be considered negligible. If fast neutrons 
occur only at one energy, the source term can be written as 

S * K Zj 



( 21 ) 
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Using Glaaston© and Edlund's (1£>) treatment the multi- 
plication factor may be simplified to 



K<» v 




s 

FUEL 



7 



( 22 ) 



by assuming that the fast fission factor, the resonance escape 
probability, and the fast neutron utilization are essentially 
unity. The fast neutron utilization term implies that the 
capture of neutrons by material other than the fissionable 
phase is negligible. Shapiro (16, p.130) states that the ra- 
diative capture cross sections for elements having mass numbers 
less than 100 and neutrons numbers less than 60 are not greater 
than £0 millibarns. If the fast neutron utilization factor 
were not unity, the treatment given by Kaplan (17, p*525) for 
thermal neutrons would still apply. 

A combination of Equations 20 to 22 yields 



7 2 J0 (?) + B 2 # (?) * 0 



(23) 



where B 2 2 ( Oj - 1)/D. 



(2k) 



For a weak absorber as described by Glasstono and Edlund 
(15>, pp.92,98), the diffusion coefficient is related to the 
transport mean free path by 
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( 25 ) 



where A^. 3 A # (1 - p^) and p * 2/3A 



( 26 ) 



Aa a first approximation the macroscopic transport and 
scattering cross-sections arc considered to be essentially 
equal except in cases where the actual values are known. The 
atomic masses of tho uranium, zirconium, and molybdenum are 
high enough to make the tern of Equation 26 negligible. As 
a second approximation the elastic scattering oross section is 
expected to be equal roughly to the total cross section in the 
middle energy region of 1 to $00 kev where neutron gamma cross 
sections ore small, Shapiro (16, p.78). The constants used in 
this analysis arc listed in Appendix D with additional refer- 



Tho solution to Equation 23 for the buckling in taras of 
spatial dimensions of a right circular cylinder as developed 
by Glasstone and Edlund (15* p.213) is 



ences. 



Substitution of Equation 2 5 into 21}. results in 



B 2 « 3 S a (7J - 1) 



(27) 




( 28 ) 



where R and H are the extrapolated radius and height of the 
cylinder in centimeters. 
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A combination of Equations 27 and 28 yields 

3*t*.<*l ' x) = + Cf ) < 29) 



Equation 29 is the approximate critical equation for an 
idealized fast homogeneous unreflected reactor in the form of 
& right circular cylinder. Appendix P shows the calculations 
that were used to determine the critical mass of uranium for 
a right circular cylinder with a physical diameter and length 
of 42 and 60 in, respectively. 

The various component weights which were based on volume 
fractions with corresponding densities are calculated in Ap- 
pendix D. 

For the flux required to produce the maximum temperature 
change, the relationship 






JLcL 



“ T 0 Jt 2 ) 



max 



G CjL * 



4EE f t (L/DJ 



(30) 



was used. The result is shown in Appendix D. 



Discussion 

The length of the heat exchanger was limited to 60 in. 
because Figure 7 showed that the maximum thrust per sq ft was 
obtained with an approximate L/D # ratio of 120. This produced 
a 0.50 in. hydraulic diameter or a 0.25 in. fuel element spacing 
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The uranium-zirconium alloy with a composition of 79*3 
per cent enriched uranium by weight was selected because of 
its 2 , 450*F melting point. A larger percentage of uranium 
would reduce this melting point which is considered a minimum 
for the high temperatures desired. The alloy composition was 
considered a constant and treated as a mixture for volume 
fractions since the actual density was unknown. 

The component volume fractions of the heat exchanger 
resulted from trial and error thickness selections for the 
uranium-zirconium and molybdenum alloys which satisfied Equa- 
tion 29 for a physical diameter and length equal to or less 
than 42 and 60 in, respectively. Since the physical length 
was fixed previously for maximum thrust per aq ft, the radius 
was considered variable. 

The nuclear reactor, when treated as a homogeneous un- 
reflected reactor according to Equation 29, was determined to 
be critical at a physical diametor of 40 in and a length of 
60 in. The large critical volume, 43*6 on ft, was due to the 
low macroscopic absorption of uranium, 0.0056 cm" 1 , which re- 
sulted indirectly from the 75 per cent volume of air. A de- 
crease in the volume fraetion of uranium would decrease both 
the total macroscopic absorption and transport cross sections 
with a resultant increase in radius if the length remained 
constant. 

The weight of enriched, 90 per cent uranium required for 
the critical volume of 43.6 cu ft was 3,622 lbs. The amount 
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of zirconium and molybdenum was 960 and 3»lj70 lbs, respective- 
ly, which resulted in a total heat exchanger weight of 8,052 
lbs# At current prices of approximately 30 dollars per gram, 
$0 million dollars would be required for the 90 per cent en- 
riched uranium. 

The flux that was necessary to generate a maximum wall 
temperature of 2200®F with a 0.0315 inch fuel element thick- 
ness was 8 x 10^ n por sq cm sec which assumed zero radial 
flux depression across the fuel element. 

Since the radial flux will follow an approximate cosine 
distribution, it is apparent that for an average maximum wall 
temperature of 2000°F which was assumed to obtain 9,500 lbs 
of thrust with a 7*7 sq ft frontal area, the fuel elements 
near the centerline would operate at temperatures exceeding 
2200®F. A proposed central control rod and reflector would 
tend to flatten the flux; however, the flux envelope resulting 
from these effects was not investigated. 

The application of the one-energy group diffusion theory 
to this analysis was a realistic approximation bocause the 
transport mean free path, 8.7 cm, was small in comparison to 
the boundary distance. The extrapolation distance, R-r, is 
not generally equal to 0.71 A^ for a fast reactor and its 
exact value can bo obtained only by cor^licated methods, Glaa- 
stone (20, p.217). However, as pointed out by Glasatone (20, 
p.2l7), there is probably no point in a fast reactor that is 
more than two scattering or transport mean free paths from a 
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boundary. This statement is supported by the large quantities 
of fuel required for the heat exchanger which has a minimum of 
11 transport mean free paths in the radial direction. 

An intermediate neutron energy reactor may offer the pos- 
sibility of reducing the fuel quantity and weight while still 
maintaining a largo volume fraction of air f or compact heat 
transfer. The slowing down density term of Equation 20 was 
neglected because air, which composed 75 per cent of the heat 
exchanger volume, has a negligible moderation effect on neu- 
trons, To continue this simplification the material for the 
fuel alloy and container was selected from elements with rela- 
tively high mass numbers. However, with the proper selection 
of lighter materials for the fuel alloy, container and reflec- 
tor, the moderation effects may increase to the extent that 
this fast reactor would evolve into an Intermediate reactor. 

It is believed that any trend toward a thermal reactor would 
decrease the fuel requirements for the configuration. 

Conclusion 

From a nuclear aspect this heat exchanger, which was com- 
posed by volume of 12,5 per cent alloyed fuel, 12,5 per cent 
molybdenum structure, and 75*0 per cent air, was or it i cal at 
a diameter and length of 1^0 and 60 inches, respectively. The 
total weight was 8,052 lbs which included 3,622 lbs of enriohed 



90 per cent uranium* Approximately 50 million dollars would 
bo required for this amount of uranium. 

The flux that was necessary to generate a maximum wall 
temperature of 22 00*F with a fuel element thickness, 0.0315 

1*5 

in* was reasonable at 10 hi per sq cm sec. 

With cost as a major consideration there is probably no 
point in a fast reactor more than several scattering or trans- 
port mean free paths from a boundary. An intermediate neutron 
energy reactor may offer the possibility of reducing the fuel 
quantity and weight for this conf iguratlon. 
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OVER-ALL FEASIBILITY CONCLUSIONS 

For an air inlot velocity of H = 0.2, a weighted radial 
average of 20Q0 e F for the maximal cylinder temperature, a 
frontal area of 7.7 sq ft, a 42 in diameter, and a 60 in 
length, the amount of heat transferred was substantial enough 
to produce 9,500 lbs of thrust with a flight Mach number of 
2.5 at 37,500 ft. 

In order to obtain an average temperature of 2000*F the 
fissionable matorial within the thin right circular cylinders 
had to reach a critical condition. This required 3,622 lbs 
of enriched 90 per cent uranium with an approximate cost of 
50 million dollars. The heterogeneous nuclear heat exchanger 
was critical with a diameter and length of I4.0 and 60 in, 
respectively and a total weight of 8,052 lbs. 

The low thrust to weight ratio cay bo compensated for by 
the negligible fuel burn up when compared to a turbojet sys- 
tem. 

Molybdenum alloy, 0.45 per cent titanium, which was pro- 
tected with a disilicide coating appeared to have the high 
temperature properties for the container material. Sandwiched 
between the molybdenum was an alloy of uraniua-sirconiusi which 
contained 79.3 per cent uranium by weight. This binary system 
was selected because of the 2450*F melting point at the above 
composition. Whether a metallurgical bond could be formed be- 
tween the molybdenum and fuel alloy was undetermined. 



36 



From the above theoretical concluaiona and Halting 
conditions, the over-all feasibility of a compact axial flow 
nuclear heat exchanger was favorable for thrust, raetal avail- 
ability, and nuclear criticality, borderline for thrust to 
weight ratio, and unfavorable for the cost of fuel. However, 
the proposed central control cor© and peripheral reflector 
should increase thrust and reduce the cost of fuel because of 
savings in critical mas. 

An intermediate neutron energy reactor may offer the pos- 
sibility of reducing the fuel requirements* With the proper 
selection of lighter materials for the fuel alloy, container, 
and reflector, the moderation effecta may Increase to the ex- 
tent that this fast reactor would evolve into an intermediate 
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HEAT TRANSFER AREA = 
it ( D| + D2)dx 



Figure 2* Nuclear heat exchanger. 
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Figure 3 



Thermal choking boundaries for a constant area 
frictionleas chamber being heated. 
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Relationship of various temperatures 
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Figure 5* Curves of stagnation tenperatures and pressures versus L/D e ratios. 
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Figure 6* The Bray ton cycle for temperature and entropy 

conditions for stations corresponding to Figure 1» 
Values for dry air are from Reference (7). 




CONSTANT WALL TEMP 

MAXIMUM WALL TEMP 



1*3 



( 2 1J/91) lSfldHl 

oo o ooooo 

oo o ooooo 

{£> in ^3" rooj — Ocr> 




(an/ni8) tt7 M- £l °M ‘Ad3VHlN3 JO 39NVH0 



Figure 7* Change of enthalpy and thrust from exit of heat exchanger to exhaust 
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APPENDIX A 



Stepwise Integration for Stagnation Temperature 
and Mach Humber as a Function of L/D Ratio 

The procedure used for integration of Equations 13 and 
14 are outlined in Reference (2) and are reproduced In part 
using the following approximations : 

1. The integral of dT^/T^ nay be approximated by 



2. The coefficients of dx/D o in Equations 13 and 14 are 
taken at ^ = (M^ + Mg)/2. 

The inlet conditions are as follows: 



For a section 2 at a distance Xg - downstream from section 
1, the following computations are carried out. 

a. Let - x-ji )/D fl * 30 

b. Estimate Mg » 0*3, 25 




H i * °* 2 

T 0,l * 68 ^ R 

Pq -y * 46*5 lb/in 2 

T 0 w * 246 0«R 



N pH * 0.664 
KrP * 0.872 
A * 1.32 

f * 0.008 



c. Calculate 2(T Q , 




^0 # l^/^0,2 * A 0,l) from Equation 



1 



13 











m* « 

•o 







w 



2 ( *0.2 * T 0.l' 

<- rp 7a 

l 0,2 ♦ A 0,l 



“Ojl w 






(1 - N ) — ^"3; M ^ 
, . K F.P } 2 a. 2 
1 + — =CT 2*~” 

1 - HP M l,2 



2Af 



X 2 “ X 1 



D 



© 



(T - 88^) 
2 __£*£ 

t O,2 + 884 



2460 , . (1 - .872) (.2) (0625) 

vr-f* '884T72- - 1 * iT;l1tbW" 



2(1.32) (.008) (30) 



2 ** 1640 



d. Calculate from Equation 15 and compare with as- 



sumed value of M.. 



- M x 2 » Vo625) * <)P ^ 2(.008) (30) * 



f 



(1 ♦ 1.4 x . 0625) (1 



.32) [ . 95 ] 



♦ 2 x 1*4 x 



» . 0625 ^| 



A Z - * 0.0500 

M 2 2 » 0.0500 + 0.0400 



K, 



0.3 



This agrees with the assumed value. 



e. If step d was not successful, repeat the previous 
steps with an improved value of 
f* After convergence is obtained, the process is re- 
peated to find the Mach nurjber at a section 3* 



APPENDIX B 



Stepwise Integration for Stagnation Temperature 

and T./D Ratio a 3 a Function of Hach Number 
' o 



The procedures used for integration of Equations II 4 . and 
17 are outlined in Reference (2) . Equation 11+ may be inte- 
grated approximately over a short Interval of Mach number by 
assuming that the influence coefficients have a constant value 
corresponding to (M^ + H^)/2 t that the magnitude of ? Q is 
constant during integration with a value of (T^ ♦ T Q 2 )/ 2 , 

and that the friction factor is constant. With these approxi- 
mations and by the substitution of Equation 17 into llj., Equa- 
tion II 4 . Integrates into 



M 2 - M 2 * 2 

* A r o,i 



- 1) 



F 

T , 



7rr-*£ * 1 

0,1 



f£ 

*"? 

( 77 2i2i) - (J£*£)-i 

1 o,i r o,i 



(17a) 



For the same inlet conditions as were given in Appendix A, 

let 

PT 2 a 0.2ip f^ 2 = 0.0576 

M x * 0,20 Mj 2 * 0,014-00 

2 * °*22 M 2 2 - Kj 2 * 0.0176 
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From Reference (2), F, r » 0.05493 and F-, * 0.003149* 

l Q 



l7a converged after three attempts to 



Equation 




1.349 or T q>2 = 1188°R 



Equation 14 gave 



X 2 ’ *1 

D 






SP 



Ln 



O.w 
0*1 " 



* 

1 O.w a 0,2 

nr» * 7-—*— 

* 0,1 * 0,1 



ST^SST 



Ln 



2.785 ■ X 

2.785 - 1.344 



SC 



16.85 



Since both Mach number and stagnation temperatures are known 
at each section, the following equations were used for comput- 
ing the remaining properties of the air: 




t q.2 (1 * J r k m i 2) 
T 0,l (1 ♦ m 2 2 ) 



* 1.34 




« 0.965 



* 1.392 
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fables 1-3 which follow have calculate values for the various 
properties of air as a function of £«/D ratios and different 
constant wall temperatures* 



fable 1, Air properties aa a function of Kach number for a 
constant wall temperature of 2660°R and a friction 
factor of ,008 
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l 0 
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p o 




•n 




•R 
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ft/aec 


lb/ft 3 


lb/in 2 


0.20 


884 
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880 
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0.1385 
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0.21* 
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30.4 


0.60 


2356 


109.52 


2182 


23.6 


1370 


0.0294 


28.6 
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Table 2. Air properties as a function of Mach number for a 
constant wall temperature of 2460*R and a friction 
factor of *008 



M 


T o 

•R 


Vo. 


T 

®R 


P 

lb/ln 2 


V 

ft/sec 


lb/ft' 


p o 

lb/ln 2 


0.20 


881* 


0 


880 


45.1 


291 


0.1385 


46.5 


0 . 214 - 


1184 


13.43 


1180 


43-4 


405 


0.0998 


44.9 


0.2Q 


1473 


29.32 


1462 


41.5 


527 


0.0769 


42.9 


0.32 


1711 


46.70 


1692 


38.9 


649 


0.0625 


40.3 


0.36 


1891 


63.98 


1863 


36.4 


768 


0.0533 


37.6 


o,4o 


2018 


79.33 


1970 


33.7 


878 


0.0465 


34.8 


0.44 


2100 


92.25 


2040 


31.2 


985 


0.0415 


32.8 


0.1*8 


2155 


102.6 


2080 


23.9 


1085 


0.0376 


31.4 


0.^2 


2192 


110.7 


2097 


26.8 


1182 


0.0346 


30.2 


0.60 


2236 


121.9 


2100 


23.2 


1363 


0.0298 


27.8 


0.70 


2258 


128.2 


2050 


19.9 


1590 


0. 0259 


25.8 
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Tabic 3. Air properties as a .function of Mach number for a 
constant wall temperature of 2060°R and a friction 
factor of .008 



M 


T o 

°R 


I,/D e 


T 

°R 


P 

lb/in 2 


p 0 

lb/in 2 


.20 


884 


0 


880 


45.1 


46.5 


.24 


1178 


17.95 


1160 


43.4 


44.6 


.28 


1437 


39.75 


1411 


41.0 


42.2 


.32 


1631 


63.55 


1592 


38.2 


39.4 


.36 


1760 


85.95 


1711 


35.2 


36.3 


.40 


1839 


104.7 


1775 


32.3 


33.3 


.44 


1886 


102.0 


1811 


29.6 


31.1 


.48 


1914 


131.9 


1830 


26.9 


28.8 


.52 


1935 


140.6 


1830 


24.3 


26.9 


.60 


1945 


145.6 


1800 


21.1 


25.3 
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APPEKPIX C 



Calculation for Outlet Stagnation Temperatures as a 
Function of L/D ft Ratios and Maximum Wall Temperatures 



With reference to Figure 2, for a heat balance on the 
fuel cylinders without axial conduction requires that 



0,w 



- T, 



Snax 

2h 



sin- 



xx 



(A) 



The longitudinal volumetric heat source term is zero when 
x *= 0, or X», which assumes zero extrapolation distance* 

Neglecting axial conduction within the fluid and any heat 
losses, a heat balance on the air in the annuli requires that 




avC D dT n 
“ p e 0 

~ilh ~ 



(B) 



Equating Equations A and D and performing the required 
integration results in 



T 0 - ■ ?!f>r (33; ) • f 1 - 0O Hr ] • <°> 

Adding Equations A and C to eliminate the stagnation tem- 
perature of the air yields 
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T 



0,w 



* 0,2 * 



Sna: 



2h 



si » 



XX 

L 



4 * 



fc 

x G c 




. (D) 



Substitution of Ec^iation 7 with A * 1 into Equation D results 
in 




(E) 



From Equation E there appears to bo a maximum which repre- 
sents the highest temperature of the fuel cylinder in the 
longitudinal direction* The point at which the maximum is at- 
tained is found by differentiating Equation E and setting the 
result equal to zero* 

- k arc tan ( - fjr (P) 

The maximum wall temperature occurring in the fuel cylin- 
der can now be obtained by Inserting Equation F into equation 
E. For purposes of computation, it is convenient to represent 
stx/L by the symbol a* 



a 5 xx/L 



(G) 




rvi 
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By substituting G into iquation F, results in 



tan a 



max 



* 

2f L 



Utilizln,. Equation H, Equation E reduces to 



T T - 2Q ™* t 

°."max °* 2 0 °p * 



(®.) [ 



1 - seca 



nax 



From Equation 17, which wa3 previously derived, 



dT, 



0 

x 0,w - l 0 



2f 



dx 



and introducing Equation A, which was simplified by 
ing h = fG0p/2, the following differential equation 
rived: 



2Q t 

jm _ max 

dT o - r sT ~ 

P © 



XX 



dx 



Integrating Equation J results in 



^Snax 11 f L 



T 0,3 T 0,2 * G C p x 






(H) 



(I) 



substitute* 
was de- 

(J) 



(K) 



However, Equation K may be combined with Equation I to elirai- 
nate t, which results in 



( 18 ) 



T n . - T n p = 2 

0,3 0,2 i - soca 

max 



The magnitude of ? n utilized the previously se- 

* w iQax 

lected wall temperatures. Since T^ ^ was t* 1 © given inlet 

conditions, a was computed for various L/D ratios result- 

max c 

ing in T n , as a function of T A and a . fable 1 of 

U,J !na ** 

Appendix C presents calculated values of T Q ^ which are shown 
in Figure 5* 



Table 4. Outlet stagnation temperatures as a function L/D 
ratios and maximum wall temperatures 



V D e 




40 


60 


80 


100 120 


140 


tm mx 




-4-910 


-3.270 


-2.450 -1.965 -1.635 


-1.405 


■•“mx 




-5.010 


-3.420 • 


-2.650 


-2.210 -1.920 


-1.728 


X -**°* M x 




6.010 


4.420 


3.650 


3.210 2.920 


2.728 






»M« 


884*R and 


^ ,w nax 


* 266 0*R 




t o,3 


Oft 


1475 


1688 


1858 


1990 2101 


2184 






T 0,2 * 


084°R and 


7 0,w 
* nax 


* 2460°R 




T 0.3 


°R 


1408 


1597 


1748 


1865 1964 


2038 






H* as 

x 0,2 


884°R »nd T n M 

;#w max 


a 206 0°R 




T o,3 


°R 


1275 


1415 


1529 


1616 1689 


1745 
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APPENDIX D 



Critical Calculations for a Homogeneous Unreflected Reactor 

The effective neutron energy was taken as 0*3 Kev based 
on flux spectra computed for existing fast reactors by Murray 
(18, p*69)* Prom the same Reference the following constants 
were used: 

D 238 : <T p * 0.18 barns <T^ * 0.0 <T, a barns 

U 2 ^ : 0" a * 1.78 barns (T « 1.50 barns <r r * 7.1}. barns 

V* 2.51 

From Hughes and Harvey (19) at 0.3 Mev, the total cross 
sections of the selected reactor materials and air were as 
follows : 



The above total cross sections were considered to be equal ap- 
proximately to the scattering or transport cross sections. 

A 90 per cent enriched uranium mixture by weight produced 



Zr : 0“ T * 9,0 barns 

Ko 1 (T^ » 9,0 barns 

>T : (TV » ij .,0 barns 

0 : <r^ = 3.0 barns 



238 Hg + 235 N5 * #9 ° 




or 



\ = 0.105 

s 




• M 

* * I • 

Ml It 





The solution to Equation 22 with the above result yielded 



f,5 



2.5 1 



Jj>jO 



cr a ,5 * *“*,8 (Kq/N^) 1.76 ♦ (.18)U05) 



2.09 



Sine® the alloyed density or the uranium-zirconium system 
was unknown, the alloy, 79*3 per cent uranium by weight, was 
treated as a mixture to obtain the Tallowing volume fractions: 

f 238 = °-°5S 

f 23S " °- 512 
f Zr « 0.1)32 

A hydraulic dianctor, D * 0.5 in. .as obtained with a 
selected L/D * 120 and L * 60 in. The following thicknesses 
were assumed: 

Uranium-zirconium alloy: 0,03125 in 
Double wall of molybdenum: 0.03125 in 
Hydraulic diameter/2: 0.25000 in 
Thickness of air space: 0.18750 in 

The approximate component volume fractions obtained from the 
above solidity factors were as follows: 

* n * 0.5925 
f 0 * 0.1575 
f 23Q * o.oo7o 

^235 55 

f Zr * 0.0540 
r Mo - 0.1250 
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The macroscopic cross sections were computed for each of 
the above volume fractions which gave for a total 

£ * 0.1163 cra" J and S * 0.0056 cm" 1 

U a 

Solution to "ruction 27 yielded 

B 2 * 21.7 x 10" 1 *- on* 2 . 

With a physical length equal to 60 in or I52.lv cm and the 
transport mean free path equal to 6.45 cm, the length includ- 
ing the extrapolation distance was 

H * h + .71 A * I52.lv + 6.0 * 156.U- cm 

The solution to Equation 29 for the critical radius in- 
cluding the extrapolation distance was 56.7 cm; the physical 
radius was 50.7 cm or 20 in. 

Since the total volume waa 12.35 x left cu cm, the follow- 
ing component weights were determined by using the previously 
computed volume fractions with corresponding densities: 



^38 


362 lbs 


u 2 35 


3260 lbs 


Zr 


960 lbs 


Mo 


3U70 lbs 


Total Weight 


8052 lbs 
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From Equation K, Appendix C, the flux that was necessary 
to generate a maximum wall temperature of 2200*? with a fuel 
element thickness, 0.031$ in, was obtained with 



Saax 



<*0-» - *0.2> 0 C T, * 

UB S r t (L/D e ) 



where 

E « 200 Mev/fiasion 

E a 0.045 cm 

Mev * 1#$2 x 10** 16 BTU 

G * 40*4 lb/ft 2 see 

C * .2671 BTU/lb •? 

P 

L/P * 120 

Tq 3 - ?q £ * 2101 - 884* from Figure 5* 

The flux obtained from the above equation was 8*25 x 10^n/ 

cm 2 sec. 
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